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Abstract 

We study the linearization of three dimensional Regge calculus around 
Euclidean metric. We provide an explicit formula for the corresponding 
quadratic form and relate it to the curl T curl operator which appears in 
the quadratic part of the Einstein-Hilbert action and also in the linear 
elasticity complex. We insert Regge metrics in a discrete version of this 
complex, equipped with densely defined and commuting interpolators. 
We show that the eigenpairs of the curl T curl operator, approximated 
using the quadratic part of the Regge action on Regge metrics, converge 
to their continuous counterparts, interpreting the computation as a non- 
conforming finite element method. 

1 Introduction 

Regge calculus [19] is a combinatorial approach to Einstein's theory of general 
relativity |61| . Space-time is represented by a simplicial complex. Given this 
simplicial complex, a finite dimensional space of metrics is defined, each metric 
being determined by a choice of edge lengths. We call such metrics Regge met- 
rics. A functional, defined on this space of metrics and mimicking the Einstein- 
Hilbert action, is provided. We call this functional the Regge action. A critical 
point of the Regge action on the space of Regge metrics is generally believed to 
be a good approximation to a true solution of Einstein's equations [H] . 

Regge calculus (RC) is quite popular in studies of quantum gravity [S^. Its 
discrete nature also makes it a natural candidate for the construction of effi- 
cient algorithms to simulate the classical field equations, a possibility expressed 
already in the last sentence of Regge's paper. In this direction we are aware of, 
in particular [45] [46] [9] [32] [31] . However, it seems that the bulk of numerical 
relativity computations are performed using other methods. One difficulty with 
simulating Einstein's equations is the gauge freedom (diffeomorphism invari- 
ance) which imposes constraints on the evolution. Hyperbolicity in this context 
is a delicate matter ^51j . Progress on the simulation of merging black holes [17] , 
seems to have been achieved in large part by judiciously choosing which partial 
differential equations to solve (in particular the gauge conditions), so that a 
number of traditional discretization philosophies, including finite difference, fi- 
nite element and collocation methods, are successfully applied today [T|[53][TT]. 
in support of the emergent field of gravitational wave astronomy |56j . 
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We hope that this paper can contribute to developing RC into a good al- 
ternative, or facilitate the integration of some of its appealing features into 
currently used methods. Its geometric "coordinate free" nature would make it a 
structure-preserving method in the sense of [33] . Thus our motivations are close 
in spirit to for instance |29j |53| . The kind of variational structure that underlies 
RC has become a governing design principle both for finite element methods 
and integration of ordinary differential equations [55] [10] , so that insights from 
these mature fields could well inspire decisive improvements in RC. 

We are not aware of any stringent convergence results for RC, except those 
of [21]. There, it is shown that for any given smooth metric, the Regge metrics 
interpolating it, have a curvature (defined by Regge calculus) which converges 
in the sense of measures, when the mesh width goes to 0, to the curvature 
(defined in the standard way) of the given smooth metric. In numerical analysis 
this would be called a consistence result. In general, consistence is only a step 
towards proving convergence of a given numerical method. We also point out 
that the convergence of RC is discussed, in less stringent terms, in the physics 
literature, e.g. |15| and references therein. 

In [35] we related the space of Regge metrics to Whitney forms [55][63j. 
As remarked in [T3], Whitney forms correspond to lowest order mixed finite 
elements [1H][13], the so-called edge and face elements, for which one has a rel- 
atively well developed convergence theory [17] [54]. More recently this analysis 
has been cast in the language of differential forms and related to Hodge theory 
[37][53][S]. We showed, in [35], that there is a natural basis for the space of 
Regge metrics expressed in terms of Whitney forms and that second order dif- 
ferential operators restricted to Courant elements (continuous piecewise affine 
functions) are in one to one correspondence with linear forms on Regge metrics, 
edge elements and Courant elements. This link integrates Regge calculus into 
the finite element framework. However we did not approach the question of 
curvature which is central to RC. 

In this paper we further develop the theory of Regge elements. We first 
insert them in a complex of spaces equipped with densely defined interpolators 
providing commuting diagrams as in finite element exterior calculus [4]. The 
differential operators of this complex are those of linear elasticity. For the 
purposes of relativity, it appears that less regularity is required of the fields 
than for continuum mechanics, so that the discrete complex we obtain differs 
from those constructed in for instance [3]: the last two spaces in our complex 
consist of matrix and vector valued measures that cannot be represented by 
integrable functions. Next we provide results concerning the Regge action. 

A priori it is not clear if RC should be considered a conforming or a non- 
conforming finite element method. Is the Regge action the restriction to Regge 
metrics of some extension by continuity of the Einstein-Hilbert action, to a 
large enough class of non-smooth metrics? One might compare with Wilson's 
lattice gauge theory discretization of the Yang-Mills equations [55], where the 
discrete action is not a simple restriction of the continuous one. Indeed, seen as 
functional, defined on the vector space of Lie-algebra valued (tensor product) 
Whitney forms, the Yang-Mills action is polynomial of order 4, whereas Wilson's 
action is transcendental, even in the case of Maxwell's equations 

Given a metric, the scalar curvature multiplied with the volume form pro- 
vided by the metric, is a certain density on space-time depending non-linearly 
on the metric (and its derivatives). For a Regge metric, which has only partial 



2 



continuity properties between simplexes, it is not clear which, if any, of the avail- 
able expressions of the curvature in terms of the metric, make sense. Partial 
derivatives of discontinuous functions can be defined as distributions or cur- 
rents, in the sense of Schwartz and de Rham, but for distributions, products are 
notoriously ill defined. In addition, if the background is only a piecewise afhne 
manifold, the associated transitions between coordinate maps are Lipschitz but 
non-smooth, so that distribution theory seems inappropriate. 

The arguments put forward by Regge to justify the definition of the Regge 
action are integral in nature. In RC there is a natural notion of parallel transport 
along paths that avoid the codimension-2 skeleton of the simplicial complex. 
Around closed loops, this parallel transport behaves as if there were a curvature, 
concentrated only on the codimension-2 skeleton, with the expression provided 
by RC in terms of deficit angles. See also the justifications provided in [50] . 
To the author, it seems desirable that this ad hoc point of view, be related 
to the contemporary mathematical theory of partial differential equations. In 
this paper we provide results concerning linearization only, but holonomics will 
nevertheless play a pivotal role. 

If we expand, as is done for instance in |60l. the Einstein-Hilbert action in 
small perturbations around Minkowski space-time, the linear term is since the 
Minkowski metric solves the Einstein equations. The first non-trivial term is 
a quadratic form which we denote by Q. The Euler-Lagrange equations corre- 
sponding to finding critical points of Q, among "all" metrics, are nothing but 
the linearized Einstein equations, which describe the propagation of infinitesimal 
gravitational ripples on flat space-time. 

It appear^ that the spatial part of Q is associated with the curl T curl oper- 
ator appearing in the linear elasticity complex ;.3|, where it encodes the Saint- 
Venant compatibility conditions. We show that this (spatial) quadratic form, 
defined a priori for smooth fields, has a natural extension to (spatial) Regge met- 
rics. Moreover we show that this extension corresponds to the quadratic part 
of the Regge action, establishing that the first non-trivial terms (the second 
variations) of the Regge action and the Einstein-Hilbert action agree. 

However, the natural Hilbert space on which the quadratic form is continu- 
ous (just) fails to contain the Regge elements! We argue that RC is a minimally 
non-conforming method. As a step towards an analysis of the convergence of 
numerical methods based on RC, we show that the eigenvalues for the curl T curl 
operator are well approximated with Regge elements. For an alternative con- 
vergence result concerning linearized RC, see [lOj . 

The theory we develop is inspired by works on the eigenvalue problem for 
Maxwell's equations, i39j|13][19] and also [23]Il]l25]. As for Maxwell's equations, 
the operator does not have a compact resolvent (due to the existence of an 
infinite dimensional kernel) , so the basic theory [7* has to be amended. Indeed it 
has been shown that in this situation, stability is not sufficient to get eigenvalue 
convergence [H]. At least two additional difficulties arise. First, linked to the 
above mentioned problem of hyperbolicity, is the fact that there are eigenvalues 
of arbitrary magnitude of both signs. One of the signs corresponds to modes 
that are excluded by the constraints in the continuous case. Due to the lack of 
sign, Cea type arguments valid for Maxwell's curl curl operator (which is positive 
semi-definite) , have to be replaced by inf sup conditions [B] |16j . Second is the 

^See the acknowledgement. 
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(limit) non-conformity of the method. Central to the argument we develop is 
an analogue for metrics of the Hodge decomposition of differential forms. 

A number of interesting related results have been published while this paper 
was under review. We mention some, that come in addition to those already 
cited. RC has been described using dual tessellations ^Tj, in a framework rem- 
iniscent of the discrete exterior calculus of [35] . In a similar vein of relating RC 
to notions of discrete mechanics, we also point out [SI], whereas [H] concerns a 
finite element point of view, as in [53] . Also of interest is [5] . Numerical methods 
based on differential forms have been studied on manifolds [55] and applied to 
general relativity in a simplified setting [52]. Regge elements have been redis- 
covered as a tool for solving equations of elasticity [22] ■ Hodge decompositions 
of tensor fields, of the type used in this paper, have been studied in Lipschitz 
domains [55] . 

Layout. The paper is organized as follows. In section [5] we study Regge ele- 
ments in finite element terms. We see what happens when we apply the Saint- 
Venant operator to them, and, based on the formula we obtain, insert them in a 
discrete elasticity complex. In section |3| we relate the linearized Regge action to 
the curl T curl operator, showing that the second variations of the Regge action 
and the Einstein Hilbert action agree. In section S) we study the discrete eigen- 
value problem for the Saint- Venant operator on Regge elements. An abstract 
framework is introduced and then applied to our case. 



2 Regge elements and linear elasticity 

Basis and degrees of freedom. We consider a space-slice S which is a 
compact flat Riemannian 3-diniensional manifold without boundary. For defi- 
niteness we actually consider: 

S = {R/hZ) X (R/hZ) X (R/hZ), (1) 

for some positive reals ^i, Z2, ^3- On S we have the Riemannian metric inherited 
from the standard Euclidean structure of R^ . 

We put V = M'^ and let M be the space of 3 x 3 real matrices. The subspace of 
M consisting of symmetric matrices is denoted §, whereas that of antisymmetric 
matrices is denoted A. Elements of V will be identified with 3x1 matrices called 
column vectors (and reals with 1x1 matrices). Thus the scalar product v ■ v' 
of two vectors v,v' S V, can also be written: 

vv' = v^v'. (2) 

We let C°°{S) denote the space of smooth real functions on S. The space 
of smooth vector fields on S can be identified with C°°(S') ® V. We regard the 
gradient of a function (at a point) to be a column vector, so that we have a 
map: 

grad ■.C°°{S)^C°° {S) (g) V. (3) 

Likewise, C°°(S') OS can be identified with the space of smooth symmetric 3x3 
matrix fields. 

We partition S into tetrahedrons by a simplicial complex Th- The set of 
/c-dimensional simplexes in 7h is denoted . As is customary, the parameter 
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h denotes the largest diameter of a simplex of Th- In 52] we will be interested 
in convergence results as /i — > 0, but until then our results concern a given h 
representing a fixed mesh. 

Regge metrics are symmetric matrix fields on S that are piecewise constant 
with respect to Th and such that for any two tetrahedrons sharing a triangle as 
a common face, the tangential-tangential component of the metric is continuous 
across the face. Thus our metrics can be degenerate - we do not impose any 
triangle inequalities, as would be necessary to define a distance from the metric. 
The continuity property imposed on Regge metrics can also be expressed by 
saying that the pullback of a metric, seen now as a bilinear form, to the interface 
between two tetrahedrons, is the same from both sides. This space of metrics 
has one degree of freedom per edge, which in the non-degenerate case can be 
taken to be the length (or length squared) of the edge, as defined by the metric. 
We proceed to give basic properties of this space, including our particular choice 
of basis and degrees of freedom. 

In [22] we related this space of metrics, which we denote by Xh, to Whitney 
forms. For each vertex a; S 7^ let A^; denote the corresponding barycentric 
coordinate map. It is nothing but the continuous piecewise affine function taking 
the value 1 at vertex x, and at other vertexes. Then the following family of 
metrics is a basis of Xh, indexed by edges e S 7^^: 

Pe = 1/2 ((gradA:,J(gradAyJ^ -I- (gradAyJ(grad A^^J^), (4) 

where the vertexes of the edge e are denoted Xe and ye- 

We define degrees of freedom as follows. For any edge e consider the 

linear form on smooth metrics: 

fie-Uh^ / iVe - XeYuiXe + s{ye - Xe)){ye - Xe)ds. (5) 

One checks that on a tetrahedron these degrees of freedom are uni-solvent on 
the constant metrics. For two edges e, e' € Tfj; we have: 

, . / if e 7^ e', , . 

Me(Pe')=|i if^^g,, (6) 

The degrees of freedom ([5]) make sense for some non-smooth metrics as well, 
in particular elements of Xh- The interpolator associated with these degrees 
of freedom is the projection 1^ onto Xh, which to a symmetric matrix field u, 
associates the unique element Uh G Xh such that: 

yeeTh ^ieiuh) ^ Peiu). (7) 

The degrees of freedom do indeed guarantee tangential-tangential continuity of 
the interpolate. 

Distributional Saint- Venant operator. Recall that the curl T curl operator 
is defined on 3 x 3 matrix fields by taking first the curl of its lines, to obtain 
a new 3x3 matrix, then transposing and then taking once again the curl of 
its lines. If one starts with a symmetric matrix field, the result is a symmetric 
matrix field. We derive an expression for curl T curl w when u € Xh- For this 
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purpose we need some expressions concerning differential operators acting on 
distributions. 

Let T be a Lipschitz domain in S, with outward pointing unit normal n G 
L°°{dT) (g) V. Let Sqt denote the Dirac surface measure on dT, defined by, for 
C°°(S'): 

(SdT,^) = / ^IdT. (8) 
JdT 

Let u be the restriction to T, of a smooth scalar or vector field on S, which 
we extend by outside T. The restriction of u to dT (seen from the inside of 
T) is denoted 7(w). By integration by parts in T, we have, as distributions on 
S: 

gradu = gradjn u — j{u)n Sqt, (9) 
cnrlu — cwcIt u + "f{u) X n SgT, (10) 
div u — divT u — 7{u) ■ n SgT- (11) 

Here, a differential operator op on the left hand side is defined in the sense of 
distributions or currents on S, whereas on the right hand side, opj^ denotes the 
corresponding operator defined classically inside T. 

Let F be a (two-dimensional) domain inside a smooth oriented hyper-surface 
of S, with piecewise smooth (one-dimensional) boundary. The oriented unit 
normal on F is denoted n and the inward pointing unit normal of dF inside 
the hyper-surface is denoted m. The Dirac surface measure on F is denoted 6p, 
while the double layer distribution is denoted S'p. The Dirac line measure on 
dF is denoted dgp- We have: 

gra,dSp — nS'p + mSgF- (12) 

This formula can be applied as follows. Let m be a smooth scalar or vector 
field on S. Recall that smooth functions can be multiplied with distributions so 
that the product uSp is well-defined (it is a distribution with support on _F). A 
Leibniz rule holds for such products. We have: 

grad(u(5F) = (gradu)(5i? + unS'p + umSgF, (13) 
cnvl(u6p) = {curlu)6p — u x nS'p — u x mSgp, (14) 
diviuSp) = (div u)Sp + u ■ n6'p + u ■ mSgp. (15) 

For any 3- vector u G V, let skeww G A be the anti-symmetric 3x3 matrix 
defined by: 

(skeww)w' ^ V X v' . (16) 

Lemma 2.1. Referring to Figure{^ consider a sector in V between two half- 
planes Fq and Fi originating from a common edge E with unit tangent t. The 
unit outward-pointing normal on the planes is denoted n. The inward pointing 
normals to the edge in the planes are denoted tuq and mi . We let Ui — ±n 
be the normal to Fi such that {mi,ni,t) is oriented. Upon relabeling we may 
suppose hq — —n and ni — n. 

The Dirac surface measure on the boundary F = Fq U Fi of the sector is 
denoted 6p and the double-layer distribution S'p. The Dirac line measure on the 
edge is denoted Sp- 
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In the sector consider a constant metric u, extended by outside of it. Then: 

curlTcurlw =(skewn)M(skewn)(5^+ (17) 
((skew rii)u(ske"w mi) — (skew no)w(skew mo))^^;. (18) 
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Figure 1: View along the edge E of a sector between half-planes Fq and Fi. 

Proof. From formula (jlOD applied to each line of u we deduce: 

curlu = u(skewn)^i?. (19) 

Transposing we get: 

Tcurlw = — (skewn)u5F, (20) 
= (skew no)M(5i?„ — (skew ni)udFi- (21) 

We now apply ([T^ to the lines of the two matrices on the right hand side and 
get the result. □ 

Remark 2.1. It is not clear from the above expression that curl T curl is a sym- 
metric matrix field. However if we denote by Rp the rotation around the vector 
t by an angle p, there is an angle 9 such that Rg sends the basis (mo, rip, i) to 
(mi,ni,i). Put: 

A{p) — {skew Rpno)u{skew Rpnio) , (22) 

so that: 

A{d) ~ A{0) = (skewni)w(skewmi) — (skewno)w(skewmo). (23) 
Taking derivatives with respect to p gives: 

A'{p) — —{skew Rpmo)u{skew RpTno) + {skew Rpno)u{skew RpUo), (24) 
which is symmetric. Therefore its integral must be symmetric. 
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Proposition 2.1. Consider now a simplicial complex Th in our domain S . 

For any edge e S , let Se denote the Dirac line measure on e and tf, the 
unit oriented tangent vector along e. Let f be a face having e as an edge. We 
let rrief be the unit vector in the face f , orthogonal to the edge e and pointing 
into the face. We let n^f be the unit vector orthogonal to the face f oriented 
such that (rUef ,nef ,t(.) is an oriented basis ofY (the vector Uef depends on e 
only for a sign). Let |u]e/ be the jump of u across face f in the order ofn^f. 

We have: 

curl T curl U = ^[ujeie^e'^e, (25) 

e 

where we sum over edges e and put: 

Me = X!'^e/Ne/'^e/, (26) 
/ 

where we sum over faces f containing the edge e. 

Proof. We use the preceding lemma. The double layer distributions cancel two 
by two, by the tangential-tangential continuity of u. We have: 

curlxcurlu = (skew n^f ) {uj e/ (skew mef)Se. (27) 

Remark that the lines of the matrix (skewne/)|it]e/ are proportional to n^f (by 
tangential-tangential continuity of u). In other words there is a vector ae/ such 
that: 

(skew Tie/ )|ule/ = CtefTllf. (28) 

Then write: 

(skew Tie/) |u|e/ (skew TTie/) — ae/Tig j (skcw TTie/ ) — —a^ftl. (29) 

This matrix has lines proportional to . It follows that the lines of the matrix 
^ J (skew rig/ ) |u] (skew ruef) must be proportional to t^. Since in addition this 
matrix is symmetric, it can be written: 

(skew Uef ) juj ef (skew nief ) ^ Setetl. (30) 

The scalar coefficient Se is determined by taking traces: 

Se = tr ^ (skew Uef ) [w] e/ (skew me/ ) , (31) 
/ 

= - tr (skew ^e/ ) |m] e/ jte nlf - n^ftl), (32) 
/ 

= - ^"-e/(skewne/)Me/ie - ^ tg (skew TT-e/ ) M e/"-e/ ' (33) 

/ / 

= X^^-^/Me/^e/, (34) 
/ 

as announced. □ 

One can compare with the approach of |36j . where tangential-tangential 
continuity is not enforced on the symmetric matrix fields. The discrete curvature 
they define is based on inter-element jumps, in an expression corresponding to 
the double layer distribution in pT)) . 
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A discrete elasticity complex. The space can be inserted in a complex 

of spaces Xj^ with < A; < 3, each equipped with a densely defined interpolator 
jk 

(0) We let X^ denote the space of continuous piecewise affine vector fields. 
It consists of the vector fields of the form: 

w = ^ Wj;Aa;, (35) 

for all choices of vectors £ V assigned to vertexes x e 7^ . We equip X^ with 
the nodal interpolator I^, which to any v G C°°(S') (8) V associates: 

/^= ^ i;(a;)A,. (36) 

(1) We put X^ = Xh, and equip it with the interpolator = Ih- For any 
element of X^, its symmetrized gradient is piecewise constant and tangential- 
tangential continuous. In other words, the deformation operator, denoted def, 
induces a map X° Xl^ . For our choice of bases of these spaces we notice that 
if e is an edge with vertexes x and y, we have for any w e V: 

He{defvXx) ^ {y - x^v. (37) 

(2) We let Xl denote the space of matrix valued edge measures of the form: 



where for each edge e, is a real number, is the unit oriented tangent vector 
to e and 5^ is the Dirac line measure on e. Proposition l2. II shows that curl T curl 
induces a map — >• X'^. The standard duality on matrix fields extends to 
a non-degenerate bilinear form on Xf^ x X^, which we denote by (■,■). Notice 
that, comparing with the definition of degrees of freedom jif. on in equation 
dni), we have for any u G X\: 



tetlSe, (38) 



{Utl5e,u) = / tluU, (39) 

J e 

= fie{u)/le, (40) 

where le is the Euclidean length of the edge e. Therefore the interpolator 
deduced from the degrees of freedom /ig satisfies for any u G C°° {S) €5 S: 

VveXf, {v,llu)^{v,u). (41) 

Just as the elements of Xf^ act as degrees of freedom for X^, the elements 
of Xl^ can be used as degrees of freedom for . The associated interpolator 
onto X"^ is determined by the property that for any u G L^(S') <E)S, I^u satisfies: 

yveXl {Ilu,v) = {u,v). (42) 

We remark that, by gU and gS]), for all u e L2(S') § and all v G C°°{S)(^S: 

{Ilu,v)^{llu,llv) = {u,llv). (43) 



9 



In other words, and are adjoints of one another. 

(3) We let Xf^ denote the space of vector vertex measures of the form: 



E 



UxSx, 



(44) 



where for each vertex a:, m^; S V is a vector and Sx is the Dirac measure attached 
to X. The elements of Xf^ act as natural degrees of freedom for X^, yielding the 
nodal interpolator of X^. The standard duality on vector fields extends 
to a non-degenerate bilinear form on X^ x denoted (•,•). We define an 
interpolator onto Xf^ by requiring, for any m G (5) (g) V: 



Vve^o {Ilu,v) = {u,v). 

As in the preceding case we remark that for all u G L^('S') 
V G C°°(S')«)V: 

{llu,v}^{llu,I?^v)^{u,llv). 

One also checks that for an edge e with unit tangent tg and vertexes x and y, 
such that tf, — [y — x)/\y — x\, we have: 



(45) 

and all smooth 
(46) 



div(tetg(5e) = te{5y - Sx). 



(47) 



This shows that the divergence operator induces a map X^ Xf^. 

That concludes the list of spaces and operators we need to form our diagram. 
The following theorem summarizes some of the above remarks, and relates the 
continuous elasticity complex [3] to a discrete one. 

Theorem 2.2. We have a commuting diagram of spaces: 



1^0 



dof 



curl T curl 



■XI 



div 



div 



■XI 



(48) 

On the lower row the linear operators are defined in the sense of distributions. 

Proof. Only commutativity remains to be proved. 

(i) For any u G C°°{S) V and any edge e G 7^^ with vertexes x, ?/, we have: 



We deduce : 



/Ze(def m) = t1{u{y) - u{x)). 



Commutation of the first square follows. 

(ii) For any u G C°°{S) ® S and any v G X^, we have: 

(/^ curl T curl u. w) = (curl T curl u, ti), 
= (curl T curl w, u), 
= (curl T curl u, /^u), 
= (curl T curl /^u,?;). 



(49) 
(50) 



(51) 
(52) 
(53) 
(54) 
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so that: 

curl T curl u = curl T curl I^u. (55) 

This proves commutation of the middle square. 

(iii) Commutation of the last square follows from (i) by duality. □ 

The discrete fields considered here have less regularity than those defined in 
[2], which are all at least square integrable, throughout the complex. Spaces 
of matrix fields adapted to a second order differential operator have also been 
considered in [58]. In their case, the differential operator is "divdiv" (extract- 
ing first the divergence, per line say, of the matrix field and then the divergence 
of the obtained vector field), for which normal- normal continuity of the ma- 
trix fields is the natural analogue of our tangential-tangential continuity. Note 
also that the two-dimensional case of the "Regge complex" was studied in [27] 
(independently of any reference to Regge calculus). 



3 Linearizing the Regge action 

The setting and notations are as in the preceding section. Our aim is to compute 
the second variation of the Regge action and relate it to the previously exhibited 
Saint- Venant operator. 

We consider first what happens around a single edge. Fix an oriented line 
(edge) in V with unit tangent t. The set-up is similar to the one of Figure [TJ 
Originating from this edge are half-planes (faces) indexed by a cyclic parameter 
/ and ordered counter-clockwise. Thus the face coming immediately after / is 
denoted / -f 1. The sector between faces / and / -I- 1 is indexed by / -I- 1/2. Let 
TO/ be the oriented unit length vector in the half-plane / which is orthogonal 
to t. Let Uf be the normal to half-plane /, so that {mf,nf,t) is an oriented 
orthonormal basis of V. 

Each sector between half-planes / and / -I- 1 is equipped with a constant 
metric M/+1/2 = M/+i/2(e) depending on a small parameter e. We suppose that 
u(e) is continuous across the half-plane / in the tangential-tangential directions. 
We suppose that Uf ^1^2(0) is the canonical Euclidean metric on V, and we will 
be particularly interested in the first derivative of Uf_^_i/2{e) with respect to e 
at e = 0, which we denote by Uj_|_i/2- 

To ease notations, the dependence upon e will be implicit in what follows. 
The derivative of a function : e H> ■0(e) at e is denoted ip'ie). Unless otherwise 
specified we only differentiate at e = and therefore write "0' — 

We want to parallel transport a vector around the edge, by a path going once 
around it and with respect to the metric u{e). In each sector parallel transport 
is trivial, but from one sector to another, say from / — 1/2 to / -I- 1/2 we denote 
by Tf the matrix of the parallel transport in the basis (m f,nf,t). It is defined 
as follows. The oriented unit normal to face / with respect to uy_i/2 is denoted 
kj , that with respect to M/-1-1/2 is denoted . The operator Tf maps the basis 
(to f,kj,t) to (to f ,k'^ ,t) . We denote by i?^ the matrix of the identity operator 
from basis {mf,nf,t) to basis {mg,ng,t). 

The holonomy from sector / — 1/2 to itself, in the basis {mf,nf,t) is defined 
to be: 

%-i/2 = i?^!^/-! • ■ • R'+% ■ ■ ■ Rj+'Tf. (56) 
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Remark that Tf is an isometry from the metric to the metric uj:j^i^2- 

Consequently Ef _i/2 is an isometry with respect to Mj_i/2- In the plane or- 
thogonal to t it is a simple rotation. The angle of this rotation does not depend 
on / and is the deficit angle associated with the edge t. We denote it by 9. 

Proposition 3.1. The derivative of the deficit angle 9 at e — is given as a 
sum of jumps: 

9' = 1/2 ^ m}(K^+i/2 - u'j_„2)nf. (57) 
/ 

Proof. Let my be the oriented unit normal to t in face /, with respect to 
or equivalently m^_|_i/2- Let Pf be the matrix of the identity from basis 
{mf,nf,t) to (m f,kj,t). Since {rhf,kj) is an orthonormal oriented basis of the 
plane orthogonal to t, with respect to the metric induced by u/_i/2 we have: 

/ cos 61 -sin 6* \ 
Ef-i/2 = Pj^\ sine cos6i P/. (58) 

V 1/ 

Differentiating this expression at e = 0, using that -P/(0) is the identity matrix, 
we get: 

/ -61' \ 
^\ 0' 0. (59) 

V 00/ 

Differentiating ([5^ at e = we obtain, since T/(0) is the identity matrix: 

E'f_y2 = Rj-,T'f_,Rff' + ■■■ + RiXR) + ---+T'f. (60) 
Let J be the canonical skew matrix: 

/O -1 0\ 
J = 1 0. (61) 
\ / 

It commutes with all i?^. We have: 

6^' = -l/2tr(Ji?}_i/2) (62) 
= -l/2 5]tr(JT;). (63) 



We determine the terms in this sum. Let be the matrix in {mg,ng,t) of 
the operator sending {nig, rig, t) to {m,g,k^ ,t). We have: 

Tg^M+{M;)-\ (64) 

Differentiating with respect to e at e = we obtain, since is the identity 
matrix at e = 0: 

T'g = {M+)' ^{M;)'. (65) 

Define reals a J , , 7^ by: 



(66) 
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Then: 

tr{JT') = {a+y^ia-y. (67) 



Differentiating the identity 



™gWg+l/2fc| = 0, (68) 

at e = yields: 

"^X+l/2«9+"^g(^g^)'-0. (69) 

Since also, by the definition (|66|: 

a+ = mlk+, (70) 

it follows that: 
Similarly we have: 



Insert these two expressions in (|67p and combine with (1631) to get the proposition. 

□ 

We consider now the general case of a simplicial complex in S. For a Regge 
metric u G Xh its Regge action TZ{u) is defined as follows. For any edge e its 
length, as determined by u, is denoted leiu) and its deficit angle Oe{u). Then, 
summing over edges, one defines |49) : 

TZ{u) =Y,&e{u)le{u). (73) 

e 

Proposition 3.2. Let Regge metrics u(e) G depend smoothly on a small 
real parameter e. We use the notations: 

n{e) = n{u{e)), 0,{e) = 0e(w(e)), le{e) = ;e(^^(e)). (74) 
We suppose that u{0) is the constant Euclidean metric. Then we have: 

7^(e) = eV4 5]ee(0)^<.(0)-Ve("'(0)) + 0{e'). (75) 

e 

Proof. We have: 

0e(O) = 0. (76) 

Consequently: 

= E«(0) + C(0)) (?e(0) + el'M) + 0{e^) (77) 

e 

= eY,0'Mm + e2^(0UO)?^(O) + l/2(?^'(O)?e(O)) + ©(e^). (78) 

e e 

It is a remarkable property of Regge calculus, proved in f?S|, that for any e: 

E^e(^K(6)=0. (79) 
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This handles the first term in (|78p . Differentiating we get in addition: 



j29':{o)ieio)+0'Mm^o- (80) 



Inserting this in the second term of (1751) we get: 

7^(e) = eV2^0^(OK(O) + O(63). (81) 

e 

Finally consider an edge e and denote its extremities by x and y. We have 
^e(O) — \y — x\ and — {y — x)/\y — x\. We can write: 

Ue) ^ {{y - xyu{^){y - x))^/' (82) 

= {\y - x\' + e{y - xru\0){y ~ x) + 0{e'))'/' (83) 

= /e(0)(l + e/2 tlu'{0)te + O(e'), (84) 

so that: 

I'M = 1/2 le{0)tlu'{0%, (85) 

= 1/2 Ze(0)-Ve(u'(0)). (86) 

This completes the proof. □ 

Now insert expressions ([57]) and PO]) in ([75]) and compare with the expres- 
sions (psj) and (PE)) for the curl T curl operator. We conclude: 

Theorem 3.3. With the above notations we have the expansion: 

7^(u(e)) = (curlTcurlw'(0),M'(0)) + ©(e^). (87) 

It can be checked that curl T curl is also the operator appearing in the spatial 
part of the second variation of the Einstein Hilbert action. It suffices to write 
curl T curl in coordinates (index notation) and compare with the expression for 
the second variation provided in, for instance, [60] . 



4 Eigenvalue approximation 

An abstract setting. In this paragraph we develop a convergence theory for 
some non-conforming eigenvalue approximations. It is a generalization of the 
point of view detailed in |25) for conforming approximations of semi-definite 
operators. Some of our arguments are closer to a discontinuous Galerkin way 
of thinking [18] but we also need to cover the case of operators with spectrum 
in both ends of the real axis, so that coercivity should be replaced by inf sup 
conditions [6j. 

Let O be a Hilbert space, with scalar product (■,•). The associated norm is 
denoted | ■ |. Let X be a dense subspace of O, which is also a Hilbert space such 
that the injection X O is continuous. The norm of X is denoted || • ||. We 
suppose that we have a symmetric continuous bilinear form a on X. We do not 
require a to be semi-definite. We are interested in the eigenvalue problem, to 
find u & X and A G M such that: 

\fveX a{u,v) ^ X{u,v). (88) 
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We define: 



W = {ueX : yveX a{u, v) = 0}, (89) 
V = {ue X : yv eW {u,v) = 0}. (90) 

These are closed subspaces of X and we have a direct sum decomposition: 

X ^V®W. (91) 

We suppose that W is closed in O and that the injection V —?' O is compact. 
Let V be the closure of y in O. Let P be the orthogonal projection in O with 
range V and kernel W. 

We suppose that the map X ~^ X* , u M- a(w, •) determines an isomorphism 
V ^V* , equivalently: 

|a(u, t;)| 

inf sup , \, I > 0. 92 

The dual of X with O as pivot space is denoted X' . Notice that we distin- 
guish it from the space of continuous linear forms on X which is denoted X* . 
For instance O is a dense subspace of X' but not of X* . The map u H> (u, •) 
defines an isomorphism X' — >■ X* . The notation Y' will be used later, with the 
same meaning, for other dense subspaces Y oi O. 

Define a bounded operator K : X' ^ X a& follows. To any u ^ X' associate 
Ku = V (iV such that: 

Vw S y a{v,w) — {u,w) . (93) 

As an operator O ^ O, K ]s compact and selfadjoint. This guarantees that 
there is an orthonormal basis of O consisting of eigenvectors of K. The non-zero 
eigenvalues are inverses of those of ((88)) . with corresponding eigenspaces. 

We now turn to the the approximation of the eigenproblem (j88p by a Galerkin 
method. 

Suppose that (X„) is a sequence of finite-dimensional subspaces of O, and 
that for each n G N we have a symmetric bilinear form a„ on X„. We solve the 
eigenvalue problems: find u G X„ and A G M such that: 

Vw G X„ an{u,v) ~ X{u,v). (94) 

Decompose as before: 

Wn = {ue X„ : VveXn a{u, v) = 0}, (95) 
Vn = {ue X„ : VveWn {u, v)=0}. (96) 

We have a direct sum decomposition: 

Xn^Vn(BWn. (97) 

Define also Kn : O — ?> O as follows. To any u ^ O associate KnU — v E Vn such 
that: 

\/weV an{v,w) = {u,w) . (98) 

Since Kn has finite rank it is compact. It is also selfadjoint. Non-zero eigenvalues 
of Kn are the inverses of those of ([M]) . with corresponding eigenspaces. 
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For the eigenpairs of Kn to converge to those of K in a natural sense, the 
foUowing should be achieved pU] : 

\\K-KJo^o^O- (99) 

Our aim now is to devise sufficient conditions for this to hold. For insights into 
some difficulties that arise in the context of discretizations of the form (j94p . 
when a has infinite dimensional kernel, see in particular [12 . 

For two real sequences (a„) and (6n), estimates of the form, there exists 
C > 0, independent of the sequences, such that for all n: 

an < Cbn, (100) 

will be written: 

a„ ^ bn or bn )?= a„. (101) 

We also use the notation: 

an « 6„, (102) 

when: 

a„ 4 bn and 6„ ^ a„. (103) 

For the applications we have in mind, Xn is not a subspace of X, but we 
can weaken the norm of X to obtain a space containing X„, without loosing 
essential compactness properties. Precise statements follow. We suppose that 
we have two more Hilbert spaces X'^ and X^ with inclusions: 

X+ d X d X- ^O, (104) 

which are continuous and have dense range. The norms of X"*" and X^ are 
denoted || • ||+ and || • ||_ respectively. We suppose that a is continuous on 
X+ xX-. We define: 

W+ = {ueX+ : VveX- a{u,v)=0}, (105) 

W-^{ueX- ■.yveX+ a{u,v)=0}, (106) 

V+ = {ueX+ : yveW+ {u,v)^0}, (107) 

V- = {ueX- -.yveW" (u,w) = 0}. (108) 

We have direct sum decompositions: 

X+ (S W+ and X' = ®W- . (109) 

We suppose that is closed in O, which yields: 

W+ = W = W-. (110) 

We suppose also that that we have the inf-sup conditions: 

inf sup '"^"'"ll >0. (Ill) 
uGF+„gv- ||m||+ 

It follows that K is bounded X^' X+ . 
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We also suppose that the injeetion — t- O is compact. By duahty the 
injection V X~' is compact, so that, by composition, K is compact as an 
operator O — >■ X+. 

Concerning our spaces Xn we suppose that they arc equipped with projec- 
tions Qn ■ O ^ Xn which are uniformly bounded O ^ O and map W to Wn- 
Moreover we suppose that 5{Wn,W) — > 0, the gap being calculated with the 
norm of O. Recall the definition: 

S{W„,W)= sup ini \u-v\/\u\. (112) 

We impose of course that for any u d O there is a sequence u„ G Xn such 
that w„ — >■ u in O. We then have Q„u u in O. 

Lemma 4.1. Consider sequences Un = Vn + Wn with Vn S Vn and Wn € Wn- 
Suppose that Un ^ u in O, and decompose u = v + 11] with v £ V and w & W. 
Then f„ — )• f and Wn ^ w in O. 

Proof. We have: 

<- |m - M„p = \v - w„p + \w - w„p - 2{v,Wn) - 2{vn,w). (113) 
Considering the right hand side we have: 

\{v,Wn)\ = \{v,Wn-{I-P)Wn)\ < \v\ \Wn\ d{Wn,W) (114) 

<\v\\Un\6{Wn,W)^Q, (115) 

and also: 

\{Vn,w)\ = \{Vn,W - QnW)\ < \Vn\ \w - QnW\ (116) 

<\Un\\w-QnW\^Q. (117) 

This gives the lemma. □ 

The spaces X„ are not necessarily subspaces of X, but always oi X~ . 
Since W is closed in X~, the following proposition states that, as n — > oo, 
^iYn-, V~) 0, the gap being calculated with the norm of X~ . 

Proposition 4.1. There is a sequence e„ ^ such that for all u &Vn: 

|u-Pu| < e„||u||_. (118) 
Proof. For any u G we have: 

U-QnPu = Qn{u-Pu)&Wn. (119) 

We write: 

\PU - QnPu\^ = \Qn{u - Pu) - {u - Pu)\^ (120) 

= \Qn{u-Pu)\^ + \u-Pu\'' -2{Qn{u-Pu),Pu), (121) 
from which it follows that: 

|u-Pu|2 < \Pu - QuPu]"" + 25{Wn,W)\Qn{u - Pu)\\Pu\. (122) 

Recall the uniform boundedness of the Qn ■ O ^ O. The first term on the right 

hand side is handled by the fact that P is compact as a map X~ — >■ O. The 
second term is handled by the gap property of Wn- □ 
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The forms a„ are required to be consistent with a in the following sense. For 
all u G there is a sequence u„ S X„ such that u„ — !■ w in O and: 

lim sup 7—r = 0. (123) 

"^°°t>GX„ \\V\\- 

Moreover we suppose that we have the following weak inf-sup condition, uniform 
in n: 

1 ^ inf sup (124) 
«ev„ |m| Hull- 
Proposition 4.2. Under the above circumstances we have discrete eigenpair 
convergence in the sense of i99\) . 

Proof, (i) Define P„ : — > X„ as follows. To any u E X^ associate P„u — 
V eVn such that: 

Vw e Ki an{v,w) — a{v,w). (125) 

It follows from (|124p that P„ is uniformly bounded X^ ^ O. To get pointwise 
convergence, pick u £ Choose u„ £ X„ converging to u in O, from the 

consistence hypothesis (jl23p . Decompose u„ = + Wn with f„ £ 14, and 
Wn G W„. By Lemma l4Tl we have u„ — > Pu in O. We also have: 

in 1^ |a„(F„M - u„, w)| |a(u, w) - a„(M„,w;)| 

|P„u - i;„| =<; sup — = sup — 0. (126) 

^,gy„ ||u;||_ yj^v„ \\w\\- 

We deduce P„u -> Pw in O. 

We remark that for u G y we have: 

yveVn aniPnKu,v) = a{Ku,v) ^ {u,v) ^ an{KnU,v), (127) 

so that: 

PnKu ^ KnU. (128) 

Moreover, as already pointed out, K is compact as an operator O — > X+. 
Combining these remarks we get: 

II^-^«I1t7^o->0- (129) 

(ii) For u eW we have: 

Vw e K a„(is:„u,w) ^ {u,v) ^ {u,v - Pv). (130) 

Using (|124p and Proposition 14. 11 we can write: 

|a„(if„u,i;)| {u,v-Pv) 
\Knu\ ^ sup r = sup ^-r (131) 

vev„ \\v\\- vev^ \\v\\- 

4en\u\. (132) 

This gives: 

WKnWw^o 0. (133) 

Since K is zero on W this concludes the proof. □ 
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= x\ 


X- 





Application to Regge calculus. We first define some function spaces tliat 
will allow us to use the preceding setting. 

We define O = L^(5') ® S, equipped with the canonical scalar product, and 
let a be the bilinear form on symmetric matrix fields, which, at least for smooth 
ones is defined by: 

a(M, w) = (curl T curl u, w). (134) 
We define for any a e [—1, 1]: 

X" = {ueO : curl T curl M e H"i+"(S') ®§}. (135) 

Fix now a 6]0, l/2[ and put: 

(136) 
(137) 
(138) 

Proposition 4.3. The required hypotheses are satisfied: 

• a defines a continuous bilinear form on X , 

• the kernel W is closed in O, 

• V is compactly embedded in O, 

• a is invertible on V x V . 
Moreover: 

• a defines a continuous bilinear form on X~^ x X~ , 

• the kernel is closed in O, 

• is compactly embedded in O, 

• the inf sup condition illl]) holds. 

Proof. We base the proof on Fourier analysis. For any ^ G V we denote by 
: V — > C the associated Fourier mode: 

F{^) : X ^ exp{i^- x). (139) 

For any a £ § we have: 

curlTcurl(aF(C)) = skew (C) a skew F(^). (140) 

For £, choose two normalized vectors •^1:^2 G V such that (^7^17^2) is an 
oriented orthogonal basis of V. We make some remarks: 

First, a G § satisfies = iff a is a linear combination of the three matrices: 

<^i(o = Cier+6a, (141) 

'T2(0=ClC2"+6?[, (142) 

<^3(0 = iiCi - 6?2'- (143) 
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We compute: 

skew(Oai(C)skew(e) = -l^MO, (144) 

skew(0a2(0skew(e)= 1^1^2(0, (145) 

skew(C)a3(C) skew(0 = -{^l^a^iO- (146) 
Second, a G § is orthogonal to these three matrices iff: 

skew(Oaskew(0 = 0. (147) 
In fact this occurs iff a can be written, for some (uniquely determined) 6 G V: 

a^be + ^b\ (148) 

Since we restricted attention to periodic boundary conditions, 1^(5") Cg) C is 
spanned by Fourier modes: 

F{k) : k G {2n/h)Z x {2n/l2)Z x {2Tr/l3)Z. (149) 

The preceding remarks then provide a Hilbertian basis of O consisting of eigen- 
vectors for curl T curl. Taking into account the characterization of Sobolev 
spaces in terms of Fourier series, all the claimed results follow. □ 

We consider a quasi-uniform sequence of simplical meshes 7h of 5, with mesh 
width /i — >• 0. Attached to Th we denote as before by Xh = the space of 
Regge metrics, and equip it with the bilinear form induced by a, defined for 
w, u G Xj^ by: 

Oft, (w, u) = (curl T curl u, w), (150) 

with respect to the canonical duality pairing (•, •) on Xf^ x X^. 

The elements of X^ are not in H^^(S') (g) § since traces on edges are not well 
defined in H^(S'), therefore Xh is not a subspace of X. One might therefore be 
reluctant to call ah the restriction of a. However, consistence of ah with a is not 
in doubt. We remark also that Wh is a subspace of W (since X^ acts as degrees 
of freedom for Xf^) so that 6{Wh, W) = 0. We will use that traces on edges are 
well defined H^+"(5) L^(e), from which it follows that Xh is in X~ . 

Proposition 4.4. The spaces Xh can be equipped with projections Qh that are 
uniformly bounded 0^0 and map W to Wh- 

Proof. We use the technique of |4j , which was based on the earlier works |S7] [33] . 

Let </) : K'^ M be smooth, non-negative with support in the unit ball and 
integral 1. Define its scaling by: 

(j)f^(x) = h-^(j){h-^x). (151) 

Let R'j^ be regularization by convolution by 4'eh ■ Remark that it commutes with 
constant coefficient differential operators. Composing with the already defined 
interpolator Ih we get an operator IhRh ■ O — > Xh- We can fix a small e such 
that for all h: 

||(id-/;,i?;,)|xJ|o->o < 1/2. (152) 

We may then put: 

Qh^ihRhlxJ-^hRh, (153) 
to get the required projection. □ 
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Proposition 4.5. We have an estimate, uniform in h: 

\fuGXh ||P«||_ =^ sup ^4^- (154) 

veXu \v\ 

Proof. We first evaluate the norm of tlie restriction operator from a tetrahedron 
T G 7^ to an edge e in norms iV-+°'{T) L'^{e). 

We let T be a reference tetrahedron of diameter 1. We suppose that the 
scaling x t-^ hx, maps the tetrahedron T to T and the edge e to e. For a 
fc-multilinear form u on T, its puUback u to T satisfies, for x gT: 

u{x) = h''u{hx). (155) 

With this definition, we remark that V commutes with the pullback by scaling 
maps. Specializing to the case A: = 2 we have the scaling estimates: 

\\u\\i^2^T) = h~'/^\\uix)h2^T), (156) 

I|Vm||l=(t) = h-^'/^WVuix^^^^r), (157) 

|«|h»(t) =/i-'/'-"|ii(a;)|H»(T). (158) 

From this we deduce: 

\Jjlut,\=h-'\jAlut,\ (159) 

^/i-i/2||u||hi+.(t). (161) 
Therefore, summing over all tetrahedrons: 

(E(/e^>*e)^)'/' =^ h-'/^uy^+.^s). (162) 

e 

By duality we get, for any family of reals G M attached to the edges e: 

\\J2vetetlS4^-^-c^S) ^ h-'/\Y.<f'^- (163) 

e e 

By scaling, for such a family [vq) wc also have: 

\\Y.v,P.\\^.^s)^h-'l\Y.^lf'"- (164) 

e e 

Pick now u G and define the family (vg) by: 

curl T curl M = ^^t;efete^e- (165) 

e 

We define v S by: 

V = ^VePe, (166) 
e 

and can now write: 

|(curlTCUrlM,Z;)| ^ I Ee ^e/e^ePe^e| 

\v\ \v\ 

e 

:^ II curlTcurlM||H-i-o(5). (169) 
This completes the proof. □ 
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Corollary 4.6. There is a lower bound uniform in h: 




> 1. 



(170) 



Proof. To obtain the uniform lower bound for the second term remark first that 
Proposition 14. 1 1 yields for Uh G Vh: 



Then apply Proposition 14.51 

The equality reflects that a map and its adjoint have the same norm. □ 

Therefore the abstract setting can be applied to prove (on toruses): 

Theorem 4.7. Linearized Regge calculus yields a convergent method to approx- 
imate the eigenpairs of the Saint- Venant operator. 
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